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ABSTRACT 

Q ■ We study the effects of a global magnetic field on viscously-rotating and 

' vertically-integrated accretion disks around compact objects using a self-similar 

treatment. We extend Akizuki & Fukue's work (2006) by discussing a general 
magnetic field with three components (?*, ip, z) in advcction-dominated accre- 
tion flows (ADAFs). We also investigate the effects of a global magnetic field 
on flows with convection. For these purposes, we first adopt a simple form of 



> 

in 

C*") | the kinematic viscosity v = ac 2 s /VLK to study magnetized ADAFs: a vertical 

in 
o . 

OO , accreted but also decreases the isothermal sound speed. Then we consider a 



and strong magnetic field, for instance, not only prevents the disk from being 



more realistic model of the kinematic viscosity v = ac s H, which makes the 
infall velocity increase but the sound speed and toroidal velocity decrease. 
We next use two methods to study magnetized flows with convection, i.e., we 
take the convectivc coefficient a c as a free parameter to discuss the effects 
of convection for simplicity. We establish the a c — a relation for magnetized 
flows using the mixing-length theory and compare this relation with the non- 
magnetized case. If a c is set as a free parameter, then \v r \ and c s increase 
for a large toroidal magnetic field, while |u r | decreases but \v v \ increases (or 
decreases) for a strong and dominated radial (or vertical) magnetic field with 
increasing a c . In addition, the magnetic field makes the a c — a relation be 
distinct from that of non-magnetized flows, and allows the p on r" 1 or p cx r~ 2 
structure for magnetized non-accreting convection-dominated accretion flows 
with a + ga c < (where g is the parameter to determine the condition of 
convective angular momentum transport). 
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1 INTRODUCTION 

Rotating accretion flows with viscosity and angular momentum transfer can be divided into 
several classes, depending on different structures and energy transfer mechanisms in the 
flows: advection-dominated accretion flows (ADAFs), advection-dominated inflow-outflows 
(ADIOs), convection-dominated accretion flows (CDAFs), neutrino-dominated accretion 
flows (NDAFs) and magnetically-dominated accretion flows (MDAFs). 

ADAFs were introduced by Ichimaru (1977) and then have been widely studied over 
thirty years. The opically-thick ADAFs with super-Eddington accretion rates were discussed 
by Abramowicz et al. (1988) in details (see also Begelman et al. 1982; Eggum et al. 1988). 
The optically-thin ADAFs with low, sub-Eddington accretion rates were discussed by Rees 
et al. (1982) and Narayan & Yi (1994, 1995a, 1995b) (see also Ichimaru 1977; Abramowicz 
et al. 1995; Gammie & Popham 1998; Popham & Gammie 1998; Wang & Zhou 1999). In 
particular, Narayan & Yi (1994) introduced self-similar solutions for ADAFs with the fixed 
ratio of the advective cooling rate to the viscous heating rate in the disk. Wang & Zhou 
(1999) solved self-similar solutions for optically-thick ADAFs. The effects of general relativity 
were considered in Gammie & Popham (1998) and Popham & Gammie (1998). 

CDAFs were presented in details in Narayan et al. (2000, hereafter NIA). They discussed 
the effects of convection on angular momentum and energy transport, and presented the 
relations between the convective coefficient a c and the classical viscosity parameter a. A 
non-accreting solution can be obtained when convection moves angular momentum inward 
and the viscosity parameter a is small. Later, a series of works have been published to discuss 
the disk structure, the MHD instability, the condition of angular momentum transport in 
CDAFs (e.g., Igumenshchev et al. 2000, 2002, 2003; Quataert & Gruzinov 2000; Narayan et 
al. 2002; Igumenshchev 2002; Lu et al. 2004; van der Swaluw et al. 2005). 

The effects of a magnetic field on the disk were also studied (see Balbus & Hawley 
1998; Kaburiki 2000; Shadmehri 2004; Meier 2005; Shadmehri & Khajenabi 2005, 2006; 
Akizuki & Fukue 2006; Ghanbari et al. 2007). Balbus & Hawley (1998) discussed the MHD 
turbulence initiated by magnetorotational instability (MRI) and its effects on the angular 
momentum transportation. Kaburaki (2000) considered an analytic model to describe the 
ADAFs with a global magnetic field and Meier (2005) considered how a turbulent and 
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magnetized disk creates a global well-ordered magnetic field, and introduced a magnetically- 
dominated flow. Shadmehri (2004) and Chanbari et al. (2007) discussed the self-similar 
structure of the magnetized ADAFs in spherical polar coordinates. Moreover, Shamehri & 
Khajenabi (2005, 2006, hereafter SK05, SK06) presented self-similar solutions of flows based 
on the vertically integrated equations. They discussed the relations between magnetic fields 
components in different directions, and mainly focused on the effects of the magnetic field 
on the disk structure. Akizuki & Fukue (2006, hereafter AF06), different from SK05 and 
SK06, emphasized an intermediate case where the magnetic force is comparable to other 
forces by assuming the physical variables in the disk only as functions of radius. However, 
they merely discussed a global toroidal magnetic field in the disk. 

In this paper, we first extend the work of AF06 by considering a general large-scale 
magnetic field in all the three components in cylindrical coordinates (r, ip, z) and then discuss 
effects of the global magnetic field on the flows with convection. We adopt the treatment 
that the flow variables are functions of the disk radius, neglect the different structure in 
the vertical direction except for the z-component momentum equation. We also discuss 
magnetized accretion flows with convection, and compare our results with those in NIA, in 
which a large-scale magnetic field is neglected. 

This paper is organized as follows: basic equations are presented in §2. We obtain self- 
similar solutions in §3 and discuss the effects of a general large-scale magnetic field on the disk 
flow. In §4 we investigate the structure and physical variables in magnetized CDAFs, and 
present the relation of the convective parameter a c and the classical viscosity parameter 
a. We adopt a more realistic form of the kinematic viscosity in §5. Our conclusions are 
presented in §6. 

2 BASIC EQUATIONS 

In this paper, we use all quantities with their usual meanings: r is the radius of the disk, 
v r and Vtp are the radial and rotation velocity, Q = v v /r is the angular velocity of the disk, 
Q,k = (GM/r 3 ) 1 / 2 is the Keplerian angular velocity, S = 2pH is the disk surface density with 
p to be the disk density and H to be the half-thickness, and c s = (p/p) 1 ^ 2 is the isothermal 
sound speed with p to be the gas pressure in the disk. 

Moreover, we consider a magnetic field in the disk with three components B r , B v and 
B z in the cylindrical coordinates (r, ip, z). We define the Alfven sound speeds c r , c v and c z in 
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three directions of the cylindrical coordinates as = B 2 /(Airp). We consider that all 
flow variables are only functions of radius r, and write basic equations, i.e., the continuity 
equation, the three components (r, tp, z) of the momentum equation and the energy equation: 

-^-(rEv r ) = 2pH, (1) 
r dr 

v r d(rv v ) _ 1 d c 2 s r3 dn\ | C v Cr | C r d / ^ \ ^ ^ 

r dr £r 2 dr y dr / r dr ^ ^ ' 

IftH - ^ (VSc.) = C - + l( ^ + C;) , (4) 
f r dc 2 s c 2 s dp _ ac 2 s r 2 ( dVt\ 2 



7 — 1 dr p dr Qk \ dr 
Here we consider the height-integrated equations using the classical a-prescription model 
with a to be the viscosity parameter, and use the Newtonian gravitational potential. In the 
mass continuity equation, we also consider the mass loss term dp/ dt. In the energy equation 
we take 7 to be the adiabatic index of the disk gas and / to measure the degree to which 
the flow is advection-dominated (NY94), and neglect the Joule heating rate. 

In AF06, a general case of viscosity r\ = pv = a Pgas(PgaB +Pmas) 1 ~ M with /1 to be a 
parameter is mentioned. If the ratio of the magnetic pressure to the gas pressure is constant 
(as assumed in the self-similar structure), the solution of the basic equations can be obtained 
with replacing a by a(l + /3) 1_At . In our paper, however, we first adopt the classical form 
v = ac 2 /Qx for simplicity in §3 and §4, in which we mainly focus on the effects of a magnetic 
field on the variables v r , v v and c s . A more realistic model requires v = ac s H with both c s 
and H as functions of the magnetic field strength. We discuss this model in §5 and compare 
it with the results in §3 and §4. 

Our equations are somewhat different from those in SK05 and SK06, since we only con- 
sider the disk variables as functions of radius r, while SK05 and SK06 discuss the magnetic 
field structure in the vertical direction. More details about the basic equations are discussed 
in Appendix A. When B r = and B z = 0, our equations switch back to the equations in 
AF06, in which only the toroidal magnetic field is considered and all the variables are taken 
to depend merely on radius r. 

In addition, we need the the three-component induction equations to measure the mag- 
netic field escaping rate: 
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B r « 0, (6) 

B z = -±(v r B z )- V -^. (8) 
dr r 



3 SELF-SIMILAR SOLUTIONS FOR ADAF 

If we assume the parameters 7 and / in the energy equation are independent of radius r, 
then we can adopt a self-similar treatment similar to NY94 and AF06, 



v r [r) = -c x a\l , (9) 



r 



v v {r)=c 2 J— 1 (10) 



r 



2 . . GM 

c s(r)=c 3 — , (11) 



r 



= ^f = 2/W*^p (12) 

where the coefficients c±, c 2 and C3 are similar to those in AF06, and /3 r , fi v and f3 z measure 
the ratio of the magnetic pressure in three directions to the gas pressure, i.e., /3 r ,tp,z — 
Pmag,r,<f,z/Pgs,s- Following AF06, we also denote the structure of the surface density E by 

£(r) = S r s . (13) 

The half-thickness of the disk still satisfies the relation H oc r and we obtain 

H(r) = H r. (14) 



Substituting self-similar relations (19l- (jl4l) to equations (Ej), (jHJ) and pj), we can obtain 
the algebraic equations of ci, c 2 and c 3 : 

" ^a 2 = c 2 2 -l- [(s - 1) + &(a - 1) + ^(s + l)]cs, (15) 
- -cic 2 a = -^"^ + 1 ) C2Ca + C3 ( s + l )\lfirP^ ( 16 ) 



+ ci. (17) 



If in the cylindrical coordinates we assume the three components of magnetic field 
B r ,<p,z > 0, then v 9 (r) can be either positive or negative, depending on the detailed magnetic 
field structure in the disk. In a particular case where B r = or B^ = 0, we can only obtain 
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the value of \ci\, but in a general case where B r B v ^ 0, we are able to determine the value 
of c 2 . 

Figures 1 and 2 show the self-similar coefficients c±, |ca| and c 3 as functions of the advec- 
tion parameter / with different (fl r , fl^, (3 Z ). We consider the disk to be radiation dominated 
with 7 = 4/3, and take s = — 1/2 (i.e., p oc r~ 3 / 2 as the common case) and the viscosity 
parameter a = 0.1, which is the widely used value. 

Figure 1 shows changes of the coefficients ci, |c 2 |, c 3 with /3 Z and P v . We neglect the radial 
magnetic fields B r , and take the parameters (j3 r , (3<p, (3 Z ) in the left three panels in Figure 1 
for (0, 1, 0), (0, 1, 1), (0, 1, 3) and (0, 1, 10), and then take the value of (3 V to be 10 in 
the right three panels. As (3 r = 0, we can only obtain \c2\ without needing to determine the 
direction of v v . The coefficients c\ and C3 increase with increasing the advection parameter 
/, but |ca| decreases monotonously as a function of / except for a strong toroidal magnetic 
field. Moreover, with the fixed ratio /3 Z , an increase of (3^ makes all the coefficients \ci\ 
become larger. Oppositely, |q| decreases with increasing f3 z . In fact, with a small radial 
magnetic field /3 r ~ 0, we can obtain an analytical solution of q from equations (IT51) - (fT7j) . 
which are similar to expressions (27)-(29) in AF06, but we should replace (1 — s)/(l + s) by 
(1 — s)(l + (3 Z )/(1 + s) and (3 by P v in those expressions instead. Also, we have c| oc /~ 1 ci, 
C3 oc ci, and ci oc ^ for large ^ and ci oc for large f3 z , all of which are consistent with 
the results in Figure 1. 

As a result, from Figure 1, we first find that a strong toroidal magnetic field leads to an 
increase of the infall velocity \v r \, rotation velocity \v v \ and isothermal sound speed c s , and 
I IV I and c s are large in the case where the disk flow is mainly advection-dominated, but the 
rotation velocity \v^\ increases with increasing / only in the case where the toroidal magnetic 
field is large enough. This conclusion is consistent with the case 1 in AF06. Second, the high 
ratio fi z decreases the value of \v r \, \v v \ and c s , which means that a strong magnetic pressure 
in the vertical direction prevents the disk matter from being accreted, and decreases the 
effect of gas pressure as accretion proceeds. 

Figure 2 shows how c±, \c2\ and c 3 change with (3 r and (3^, where we neglect the vertical 
magnetic B z . For a small value of /3 r , the coefficients Ci, \o^ \ and C3 also increase with increas- 
ing (3^. However, a change of Cj is not obvious for a large value of {3 r . We are able to calculate 
the limiting value of q in the extreme case where j3 r is large enough and p r P^ 7^ using an 



analytical method. From equations (fT5l)-(fT71). we can obtain c\ = 2/[e" + y(e") 2 + 2a 2 ] and 
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c\ = e"ci for large /3 r , where e" = ^ {(7 — l)" 1 + s — 1}. If e" ^> a, we have C\ ~ (e") -1 ^ / 
and \c2\ ~ 1, which means that the infall velocity \v T \ increases with advection parameter 
/ linearly, and the radial velocity \v^\ is nearly the Keplerian velocity, no matter whether 
the disk is efficiently cooled or not. Also, for a large value of j3 T , equation ( fill becomes 



— (cic 2 a)/2 ~ c 3 (s + l)J/3 r f3<p. Since ci, c 3 > in the accretion disk, we obtain c 2 < 0, which 
means that the direction of rotation in the disk is opposite to the toroidal magnetic field 
Bp. Actually, in the case where (3 r is sufficient large and (3 r f3 Lp ^ 0, the angular momentum 
transported due to the magnetic field stress is dominated over that due to the viscosity, and 
balances with the advection angular momentum. As we take B r ^^ z > 0, from equation (j3j), 
we obtain that the large angular momentum due to the magnetic field stress makes the value 
of the advection angular momentum (Mv^r, where M is the mass accretion rate) increase 
in the disk, which requires v v < in the self-similar structure^. 

From the mass-continuity equation (TjQ) and the induction equations (J6])- (JHJ) as well as the 
solved coefficients q, we can solve the self-similar structure of the mass loss and magnetic 
field escaping rate with forms of p = p r s ~ 5 ^ 2 and B ri(PjZ = B romhZQ A s ~ 5 ^ 2 where p satisfies 



Po = -{ S+ 2) 2H • (18) 
As mentioned in AF06, when s = —1/2, i.e., S oc r -1 / 2 or p oc r~ 3//2 , there is no wind in the 
disk, and thus we can use the formula — 2irrv r T l = M to determine the surface density E. In 
addition, in the region of the disk where is adiabatic with p oc r" 1 ^ 7 " 1 ), p oc r -7 ^ 7-1 -', v r oc 
r (3-27)/(7-i) ([ e ^ s — [j_2)/( , y— 1)), we obtain the self-similar solution of c\ol = V2,C2 = 
and C3 = 0, which describe the Bondi accretion. However, if the disk region satisfies the 
entropy- conservation condition with / = 0, we can still obtain an accretion-disk solution 
beyond the self-similar treatment. For a CDAF with p oc r" 1 / 2 (NIA), a steady disk without 
wind requires C\ = or v r = 0. 
B r0t[ p 0tZo satisfy 

B ro « 0, (19) 



1 In some previous works (e.g., Wang 95, Lai 98, SK05 and SK06), the rotation velocity is taken to be positive and the toroidal 
magnetic field B v to be negative. The advection transports angular momentum inward, while the magnetic stress transports 
angular momentum outward instead. This previous result is consistent with our result here if we change the cylindrical coordinate 
used above from (r, (p, z) to (r, — tp, z). In this paper it is convenient for us to take B T)V ,z > and to obtain s series of self-similar 
solutions about magnetized flows in many different cases. 
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where ifo i n the expression ()14p can be obtained from the hydrostatic equilibrium equation 
(HI), that is, 



H = - 
u 2 



(s - l)c 3V /A^ + ^/c§(s - l) 2 (3 r (3 z + 4(1 + J v + /3 r )c 3 



(22) 



and thus we obtain the half-thickness of the disk H = H$c s / '{tJc^Ik)- We will discuss the 
effects of the magnetic field on H later in §5. 



4 SELF-SIMILAR SOLUTIONS FOR CDAF 

In CADFs, both advection and convection play contributions to the angular momentum 
and energy transportation. We propose a CDAF model in a global magnetic field in order to 
compare it with non-globally-magnetized CDAFs. We follow the idea of NIA in this section 
and consider the effect of a magnetic fielcR The viscosity angular momentum flux is 

-2 



J "— a w K pr3 ^ (23) 

and the convection angular momentum flux can be written as 

j c = _« c ^L pr 3(i+ 9 )/2^_ f nr m- 9 )/A ; (24) 



' Qk dr 

where a c is the dimensionless coefficient to measure the strength of convective diffusion, 
g is the parameter to determine the condition of convective angular momentum transport. 
Convection transports angular momentum inward (or outward) for g < (or > 0). 
The energy equation of a CADF is 

^ + ^(^) = ^ = / i5± ^(f)'. <*) 

where the convective energy flux F c is 

where we still consider the general energy equation without vertical integration as in NIA, 
and still neglect the Joule heating rate. 

2 We adopt the (a, a c )-prescription following NIA. The MHD simulations beyond this prescription can be seen in Igumenshchcv 
et al. (2002, 2003), Hawley & Balbus (2002) and so on. Moreover, Quataert & Gruzinov (2000) also develop an analytical model 
for CDAFs. 
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Using the angular momentum equation, the energy equation of the CDAF and the self- 
similar structure (l9l- lfH|) . we can obtain the relevant algebraic equations of self-similar 
structure of the CDAF, 



-cic 2 a = —^{ol + ga c )(s + l)c 2 c 3 + c 3 (s + 1)^/3^, (27) 

J ~ D (t^T + s ~ X ) ° 3ac + (t^T + s ~ x ) Cl0t = ^ a + gac " > ~4 C ^ ( 28 ) 

The radial momentum equation is still the same as that in ADAF. Combining equation ([TBI) . 
( ]27|) and (1281) . we can finally solve the coefficients c\, c 2 , and C3 in the case of CDAF and 
compare them with those in ADAF. The dimensionless coefficient a c can be calculated using 
the mixing length theory, and we adopt equation (15) in NIA, who describes the relation of 
a c with s, C3 and 7 (in NIA, they used the symbols a and Co, where a = 1 — s, and c\ = C3 
in our paper). 

To simplify the problem, we first prefer using a simpler treatment with a fixed a c to 
discussing the solutions of (1T5|) . ( 1271) and ( l28l) . i.e., we take a c as a free parameter rather 
than a calculated variable, since a c does not dramatically change in many cases. Then we 
can adopt a similar treatment as in §3 to solve equations (fl5|) . (1271) and (|28|) . Similarly as 
in §3, we first use the analytical method to discuss some particular cases. 

When {3 r f3 v ~ (which implies that the radial or toroidal magnetic fields are weak), 
we can obtain an analytical solution similar to that in §3 (see Appendix B for more details 
about the calculation). We discuss two cases. One is that the toroidal magnetic field is 
dominated and the radial magnetic field is weak (B r w 0). Then we can have an approximate 
solution, 

oict ~ - 2 _f" - y (29) 
6{a + ga c ) 



2 CiOi 

c 2 " 



a + ga c 



3 (a + ga c )(s + 1) 



(30) 



3 9(a + ga c y(s + l)' 1 ; 

where £" = ^j(s — ^("Zi + s — 1). From these equations, we know the coefficients c\ and c 3 
increase with increasing the convective parameter a c for g < (i.e., convection transports 
angular momentum inward), but we cannot obtain the relation between a c and \c%\ unless 
the values of s, g and 7 are given in detail. The other case is that the vertical magnetic field 
B z is dominated. In this case, we obtain 
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3(a + ga c )(l + s) 



3(1 



s)f3 z 



i 



\e\ 



3(a + ga c )(s + 1) 



c 3 



(32) 
(33) 

a-)A- (34) 

We find that the coefficients C\ and | C2 j decrease with increasing a c for g < while the value 
of C3 is more or less the same for a fixed /3 2 . 

Another analytical solution can be obtained when (3 T is large and (3^ 7^ 0, and we have 
the relations 



Ci oc 



a 



a + ga c 



\ 



a 



a + ga c 



|c 2 | oc 

C 3 OC -CiC 2 



2(a + ga c 



-1/2 



2a 2 



(35) 

(36) 
(37) 



From formulae fl35l) and (136]) . we again find changes of c\ and |c 2 | with a c , which is similar 
to the former case. From equation (|37|) and Ci^ > 0, we get c 2 < 0, which has also been 
obtained in §3. 

Figure 3 shows some examples of the effect of the convection parameter a c on the three 
coefficients q. In order to see the results clearly, we take a = 1 and change the value of a c 
from to 0.9 with several sets of magnetic field parameters (@ r , (3^, (3 Z )= (0,3,0), (3,3,0) 
and (0,0,3). Also we set 7 = 4/3, s = —1/2 and g = —1/3. The basic results in Figure 
3 are consistent with the above discussion using the analytical method. In particular, we 
notice that the three coefficients do not change dramatically in the case of (3 r ~ (3^, since 
the magnetic field gives a contribution to the angular momentum rather than the viscosity, 
and reduces the effect of convection on the disk. 

Next we want to obtain the relation between a c and a following NIA, i.e., we consider 
a c to be the variable as a function of s, 7 and C3. Using the treatment in NIA based on the 
mixing length theory and equations (fT5|) . (l27j) and (1281) . we can establish the a c -a relation. 
NIA discussed such a relation with g = 1 and g = —1/3, and find that the solution with 
s = —1/2 is available only for a greater than a certain critical a cr i t when the isothermal 
sound speed reaches its maximum value, and the value of a c decreases monotonously as a 
increases. However, our results are quite different from those in NIA for two reasons. First, 
we keep the term v r dv r /dr in the radial momentum equation, while in NIA this term is 
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neglected. As a result, in many cases, the sound speed to determine the actual critical a for 
available solutions does not reach exactly its maximum value. Second and more importantly, 
we consider the effect of the large-scale magnetic field on the disk. 
From equation (fl5|) . we obtain 

X -c\c? + [(1 - s)(l + /3 Z ) - /y 1 + s)]c 3 - 1< 0. (38) 

If the radial magnetic field is weak, then we have C3 < [(1 — s)j3 z }~ 1 for a large vertical 
magnetic field and C3 < 2/3 v /[9(a + ga; c ) 2 (l + s)] for a large toroidal magnetic field. On the 
other hand, from NIA, we have C3 > 7/ [(2 — s)(2 + 57 — s)] for the convective process to be 
available. Therefore, the structure of flows with convection cannot be maintained for a large 
vertical magnetic field. Moreover, if the term f3 r P v is large, we still obtain a small value of 
C3 ~ ci|c2|a / '(2 J ' PrP^) with a small value of a c , which is almost independent of the variation 
of a. 

Figure 4 shows examples of the a c — a relation with different magnetic field structures. 
We take 7 = 1.4 and s = —1/2. The left panel shows the a c — a relation with different 
values of f3 v . When the magnetic field is small (/3 ¥ ,=0 and 1 in this panel), a c decreases 
with increasing a, and a has its critical (minimum) value for the solution to be available. 
These results are basically consistent with those in NIA. However, when (3 V becomes large, 
a c increases as the viscosity parameter a increases, and the critical value of a becomes 
extremely small or even disappears. The right panel of Figure 4 shows the a c — a relation 
with different values of (3 Z . When f3 z becomes large, the critical value of a also disappears, 
but a has its maximum value. This result is quite different from NIA, who found that only 
the minimum value of a exists and becomes important. 

From the above discussion, we conclude that a strong vertical magnetic field or large f3 r f3 v 
prevents the convective process in flows, while a moderate vertical magnetic field is available 
for small a. A strong toroidal magnetic field with weak radial field makes the convective 
process become important even for large a in flows. 

In NIA, a self-similar convection-driven non-accreting solution with s — 1/2 (i.e. S oc 
r 1//2 ) was given for a + ga c = when a is smaller than the critical value a cr i t . However, 
the relation a + ga c = cannot be satisfied if (3 r Ptp 7^ for magnetized CDAFs. In fact, we 
are still able to get a self-similar structure for magnetized CDAFs when the a c — a relation 
mentioned is no longer satisfied (i.e., inequality (1581 is not satisfied). For /3 r /3 v 7^ 0, the zero 



12 Dong Zhang and Z. G. Dai 

infall velocity (i.e. c\ = 0) requires s = — 1 (p oc r -2 ) from equation ( 1271) . and a as a function 

of c 3 , 

a = a c \\g\ -|— I < a c \g\, (39) 

with the maximum value of a c to be 

(l + ft) | 9 - (2ft + 5) 7 .... 
a -'" 2 = ^7T^ 2 7 (1 + ft) ' (40) 

Furthermore, if we turn the radial momentum equation from its vertical integration to 

its general from, we can still have the relation (1391) but s = (p oc r _1 ), and the maximum 

value of a c to be 

_ (1+&) | ^(l + &) 7 

ac ' crit2 " ~wr^ 7(1+ ^ ' ( } 

Such a structure of p oc r^ 1 was also obtained by Igumenshchev et al. (2003), who explained 
the structure as a result of vertical leakage of convective energy flux from the disk. In our 
model, however, we show that this structure is due to the inefficient angular momentum 
transfer by viscosity and the zero Lorentz force in the (^-direction. 

As a result, we obtained a self-similar solution for magnetized CDAFs with c\ = 0, 
s = — 1 or s = (for the general form) and a + ga c < 0. This solution is adopted when the 
normal self-similar solutions mentioned above for convective flows cannot be satisfied. 



5 A MORE REALISTIC FORM OF KINEMATIC VISCOSITY 

In the above sections §3 and §4, we assume the kinematic viscosity v = ac^/Qx and take the 
viscosity parameter a as a constant in our discussion for simplicity. A more realistic model 
based on the physical meaning of the viscosity parameter is v = ac s H with H ^ c s /Qk in 
the magnetized disk. In this section we consider the effect of different forms of kinematic 
viscosity v. In order to compare with the results in §3 and §4, we replace a in the last 
two sections by a' and take v = ac s H = a'c^/Qx in this section. Also, we still adopt the 
definition of c\ using equation (jnj). From formula ( 1221) . we are able to obtain 



a' = a 



. 2 

S 



PrP z C 3 + (1 +p r + 



I 1 / 2 ,1 



[PrM 11 * } ■ ( 42 ) 



When (3 TyZ = 0, we have a 1 = aJl + (3^ and equation ( 1421) switches back to the case of 



fx = 1/2 in AF06. For large (3 r or j3 v and small (3 Z , we have a' ~ ay 1 + (3 r + (3^ and a' ^> a. 
For large (3 Z , we obtain a' ~ a(l + f3 r + (3^)/ (1 — s)V ' [3 r [3 z cz and a' <^ a. This result can be 
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explained as being due to the fact that a large toroidal or radial magnetic field makes the 
disk half-thickness H become large and increase the kinematic viscosity (since v oc H) , but 
a large vertical field reduces the height H and decreases the kinematic viscosity. 

Figure 5 shows the effect of a modified kinematic viscosity on the three coefficients ci, 
| C2 j and C3 in ADAFs. A more realistic expression of v increases the infall velocity, but 
decreases the radial velocity and the isothermal sound speed. However, a difference between 
these two cases of kinematic viscosity is obvious for a large toroidal magnetic field rather 
than a large vertical field. In fact, if the toroidal magnetic field B v is strong and dominated 
in (B r , 13^, B z ), we can adopt a similar solution of (49)-(51) in AF06 for fi = 1/2, and 
find that c\ increases but and c 3 reaches their limiting values with increasing (3^,. If the 
radial magnetic field B r is strong and dominated, we can obtain c\ ~ const, | C2] oc (3~ l / A 
and C3 oc /J^ 1 / 2 , which are different from §3 in which \c%\ ~ 1 for large (3 r . Furthermore, 
if j3 z is large enough, we have the limiting value c x ~ 3(1 + (3 r + /5 </3 )(l — s) -3 / 2 /?" 1 /?" 1 / 2 , 
c 3 ~ (1 — s)^ 1 /?" 1 and c\ = 3e"(s + l)c 3 , and the values of |c 2 | and c 3 are more or less the 
same, no matter what the form of kinematic viscosity is. 

For flows with convection, it is convenient for us to adopt the general definition of a c from 
NIA, which measures a degree of convection in the flows. We find that the conclusions in §4 
are not basically changed if we replace a in §4 by a'. The a' — a c relation can be turned back 
to the a — a c relation using equation (1421) . However, there is no dramatic change between 
these two relations except for extremely strong magnetic fields. 

6 CONCLUSIONS 

In this paper we have studied the effects of a global magnetic field on viscously-rotating and 
vertically-integrated accretion disks around compact objects using a self-similar treatment. 
Our conclusions are listed as follows: 

(1) We have extended Akizuki and Fukue's self-similar solutions (2006) by considering a 
three-component magnetic field B r , B^, and B z in ADAFs. If we set the kinematic viscosity 
v = ac^/flx as its classical form, then with the flow to be advection-dominated, the infall 
velocity \v r \ and the isothermal sound speed c s increase, and even the radial velocity 
can exceed the Keplerian velocity with a strong toroidal magnetic field. The strong magnetic 
field in the vertical direction prevents the disk from being accreted, and decreases the effect 
of the gas pressure. For a large radial magnetic field, v r , and c s can reach their limiting 
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values, and the direction of radial velocity is actually negative, since the angular momentum 
transfer due to the magnetic field stress in this case is dominated over that due to the 
viscosity in the disk, and makes the value of advection angular momentum increase inward. 

(2) If the convective coefficient a c in flows is set as a free parameter, \v r \ and c s increase 
with increasing a c for large and weak B r . Also, \v r \ becomes smaller and \v v \ becomes 
larger (or smaller) with increasing a c for a strong and dominated radial (or vertical) magnetic 
field. 

(3) The a c — a relation in the magnetized disk is different from that in the non-magnetized 
disk. For large B v and weak B r , a c increases with increasing a, the critical value a cr i t to 
determine different cases of the a c — a relation disappears, and E oc r~ 1//2 can be satisfied 
for any value of a. A moderate vertical magnetic field is available for small a. The large B z 
or B r B v , on the other hand, prevents the convective process in flows. 

(4) The self-similar convection envelope solution in NIA should be replaced by C\ = 0, 
a + ga c < and s = — 1 (p oc r~ 2 ) for the vertical integration form of angular equations and 
s = (p oc r _1 ) for the general form in magnetized CDAFs. This solution can be adopted in 
the region that does not satisfy the normal self-similar solutions for flows with convection 
and a c < a CiCrit 2. 

(5) The magnetic field increases the disk height H for large B r and B v , but decreases 
it for large B z in the magnetized disk. A more realistic model of the kinematic viscosity 
v = ac s H makes the infall velocity in ADAFs increase and the sound speed and toroidal 
velocity decrease compared with the simple case when the form v = ac^/Qx is assumed. 
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APPENDIX A: 

The momentum equation of accretion flows can be written as (Frank et al. 2002) 

(v ■ V)v = — VP - V$ + Q 2 r + (V • a) + — j x B, (Al) 
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where a is the viscosity stress tensor, j x B/(pc) is the density Lorentz force. Also, the 
Ampere's law and the induction equation (Faraday's law) are 

j = ^(VxB). (A2) 
<9B 

_ = V x (v x B) + Vm V 2 B. (A3) 

where rj m = c 2 /(4iro- e ) is the magnetic diffusivity and o~ e is the electrical conductivity. For 
simplicity, we consider the extreme case that a e — ► oo and r\ m m , and then neglect the 
second term in the right side of the induction equation (1A3I) . Combining equations (1A1I) and 
(1A2[) . we can obtain the three components of the momentum equation. In particular, the 
three components of the Lorentz force in the cylindrical coordinates are 

— x B)„ = -B^ + B r ^ + B z ^, (A5) 
c r or dz 

fa x B), = + BJ) + B^. (A6) 

Based on the consideration that all flow variables including the magnetic field are mainly 
functions of radius r, we can conclude v z = and d/dz = 0. Or a more realistic consideration 
requires d/dz ~ (H/r)d/dr d/dr. Also, we take d/dip = for the axisymmetric disk. We 
rewrite the Lorentz force using the Alfven sound speed as 

^ jxB >--^« + ^-7> (A7 » 

1 led 

— (j x B)^ = -c^c r + — (v^oc^), (A8) 
pc r y/por 

l(j x B). = ^|(v^J, (A9) 

These expressions are different from SK05 and SK06, who considered the magnetic field 
structure as a function of both radius r and height z: B r (r,z) = z{B t )h/ H, B^r, z) = 
z(B v )h/H with H to be the half-thickness of the disk, and B z (r,z) = B z {r). However, 
in this paper we take the magnetic field to be homogeneous in the vertical direction and 
neglect the term of d/dz as mentioned above except for the z-component equation, in which 
we take the total pressure as Ptot = -Pgas + (B 2 + B 2 )/8ir in the vertical direction, and 
adopt dPtot/dz ~ —P tot /H to estimate the value of H. In §2, we use the height- integration 
equations. 
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APPENDIX B: 



Expressions (29)- (37) in §4 can be derived as follows: 

When j3 T = or j3 v = 0, we obtain equations for the three coefficients q in CDAFs as 
1 



2 2 



[(S-1)(1 +/ 3 Z ) + (1 + S )/^]C3, 



(Bl) 
(B2) 
(B3) 

a + ga c a + ga c 

with e" = iffi^zi 4- s — 1) and £" = ^j(s — ^("Zi + s — 1). Then we can write the equation 
for c\ as 



Ci« = 3 (a + ga c )(s + l)c 3 , 



c 2 — e : , ~ +s 



^c?a 2 + da I e" + ^' 



1 

3(s + l)(a + ga c ) 3 



1 



-1 = 0. 



When (3 Z is large, the above equation can be simplified as 



loo Cia ( 1 — s 
2 C ^ + "3 



1 + sj a + ga c 



0. 



and we obtain 

3 (a + gtt c )(l + s ) 



cia 



(B4) 



(B5) 



(B6) 



(1 - s)(3 z 

Similarly, we can get the solution for large (3^ and small (3 r . 

On the other hand, if the radial magnetic field is strong and dominated and (3 V ^ 0, then 
equation ( 1B2I) should be replaced by 



- ~cic 2 a = (s + lfayj r v , 
and we obtain an equation in the extreme case, 



2^,2 



—c x a 



ac\ 



2 1 a + ga c 

and get 

ac\ = 2 (a + ga c ) 



-1 = 0, 



-112 



2(a + ga c 



(B7) 



(B8) 



(B9) 
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Figure 1. The self-similar coefficients ci, |c2 1 , and C3 as functions of the advection parameter / for different sets of parameters 
/3r, ji<p and f3 z . We take c?=0.1, 7 = 4/3 and s = —1/2. The left three panels correspond to (/3 n (9^, /3 Z )= (0, 1, 0) (solid lines), 
(0,1,1) (dashed lines), (0,1,3) (dotted lines), and (0,1,10) (dash-dotted lines). The right three panels correspond to (f3 r , (3tp 
Pz)= (0,10,0) (solid lines), (0,10,1) (dashed lines), (0,10,3) (dotted lines), and (0,10,10) (dash-dotted lines). 



Self-similar flows with magnetic fields 19 




Figure 2. The self-similar coefficients ci, |ca | , and C3 as functions of the advection parameter / for different sets of parameters 
j3 r , 0tp and /3 Z with »=0.1, 7 = 4/3 and s = —1/2. The left three panels correspond to (/3 r , f3 v Z )= (0,1,0) (solid lines), 
(0.1, 1, 0) (dashed lines), (1, 1, 0) (dotted lines), and (10, 1, 0) (dash-dotted lines). The right three panels correspond to (/3 r , f) v 
/3 Z )=(0, 2.5,0) (solid lines), (0.1,2.5,0) (dashed lines), (1,2.5,0) (dotted lines), and (10,2.5,0) (dash-dotted lines). 
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Figure 3. The coefficients ci, | C2 1 , and C3 as functions of / with different sets of parameters a c and (f3 r , f} v (3 Z ). We take a=l, 
7 = 4/3, s = —1/2 and g = —1/3. The left three panels correspond to (f3 r , P v /9 Z )= (0,3,0), the middle three panels to (f3 r , 
Pip As)= (3,3,0), and the right panels to (/3 r , fi v f3 z )= (0,0,3). Different lines refer to a c = (solid lines), a c = 0.3 (dashed 
lines) and a c = 0.9 (dotted lines). 
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Figure 4. The convective coefficient ct c as a function of viscosity parameter a with s = —1/2, 7=1.4, / = 1 and different sets 
of parameters (fir, v ,j3z)- (a) Left panel: (/3 r , /3^,/3 z )= (0,0,0) (solid line), (0,1,0) (dashed line), (0,3,0) (dotted line) and 
(0,5,0) (dash-dotted line); (b) Right panel: (/3 r , P v ,Pz)= (0,0,0.1) (solid line), (0,0,0.5) (dashed line), (0,0,0.7) (dotted line) 
and (0, 0, 0.8) (dash-dotted line) 
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Figure 5. Comparison between two forms of the kinematic viscosity u, while the results of u = ac^/flx are shown by thin 
lines, and those of v = ctc s H are shown by thick lines. We adopt s = —1/2, 7 = 4/3 and a = 1. (a) Left panels: (/3 r , (3 V 
f3z)= (2,0,0) (solid lines) and (2,2.5,0) (dashed lines); (b) Right panels: (/3 r , f3 v /3 Z )= (2,0,1) (solid lines) and (2,0,5) (dashed 
lines). 



